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We construct optimized implementations of the CNOT and other universal two-qubit gates that, 
unlike many of the previously proposed protocols, are carried out in a single step. The new protocols 
require tunable inter-qubit couplings but, in return, show a significant improvements in the quality 
of gate operations. Our optimization procedure can be further extended to the combinations of 
elementary two-qubit as well as irreducible many-qubit gates. 



According to one of the central results obtained in quantum information theory, an arbitrarily complex quantum 

■ protocol can be decomposed into a sequence of single-qubit rotations and two-qubit gates [1]. However, despite its 
. providing a convenient means of designing logical circuits, in practice such a decomposition may not necessarily achieve 
' the shortest possible times of operation and, consequently, the lowest possible decoherence rates. In any practical 

O ' implementations, however, the latter are crucially important, for any realistic qubit system would always suffer a 
^ ' detrimental effect of its dissipative environment. 

\—A . Recently, there have been various attempts to improve the performance of universal quantum gates by searching for 
' their optimal implementations among the entire class of two-qubit Hamiltonians with the most general time dependent 
I , coefficients. However, a typical outcome of such a tour-de-force variational search [2] tends to be a complicated 
sequence of highly irregular pulses whose physical content often remains largely obscure. 

In the search of a more sophisticated analytical approach, a number of authors invoked optimal control theory with 
the goal of implementing a desired unitary transformation independently of the initial state. The resulting complex 
QQ , system of nonlinear integral-differential equations can be solved numerically with the help of the Krotov or similar 
' iterative algorithms [3]. 

Conceivably, a significantly simpler alternative to the above approaches would be a straightforward implementation 
of a given unitary transformation in the smallest possible number of steps, during each of which the Hamiltonian 
remains constant. A well known example of this kind is provided by the two-qubit SWAP gate which can be readily 
implemented (up to a global phase) with the use of the spin-rotationally invariant Heisenberg inter-qubit coupling 
that remains constant during the gate operation. 

In this paper, we construct one-step implementations of some widely used universal gates. In contrast to the previous 
works where a constant decoherence rate was assumed and, therefore, the overall loss of coherence accumulated during 

■ a gate operation would be evaluated solely on the basis of its total duration, we quantify the adverse effect of the 
^ , environment by actually solving the corresponding master equation for the density matrix of the coupled qubits. In 

this way, we account for the fact that the decoherence rates generally depend upon (and vary with) the adjustable 
^ , parameters of the Hamiltonian. 

The problem of implementing a given unitary transformation in the course of quantum mechanical evolution of a 
generic two-qubit system can be formulated as the condition of a minimum deviation 
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|X-Texp(-i / H(x^(t'))dt')\\^min (1) 







between the time-ordered evolution operator governed by the Hamiltonian H and the target unitary transformation 
X. Here Xi{t) represent the control parameters whose physically attainable values are generally bounded (|xi(t)| < a^), 



and the deviation measure is given by the Frobenius trace norm ||y|| — \J Tr\Y'^Y]. 
In the case of interest, the Hamiltonian 

H^Y. ^^(B,{t) + h,{t)) + J2 Ja{t)(Jia^2 (2) 

i— 1,2 a—x.y,z 

operates on a direct product of two (pseudo)spin Hilbert spaces. 

The tunable local fields Bi = (Ai, 0, e^) encode the tunneling amplitudes Ai_2(i) and biases ei,2(i) of the individual 
qubits, while their random (and uncorrelated, < hi(t)h2(t') >= 0) counterparts hi = (0,0, /i^) represent the noisy 
environment in the worst-case scenario of two independent dissipative reservoirs. As compared to the opposite limit 
of two strongly correlated reservoirs (collective decoherence), in this case no decoherence-free subspace can possibly 
exist. 
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In the standard basis where erf 2 are diagonal, the noiseless part of Eq.(2) takes the form 

(3) 



( -'^z + ei + £2 A2 Ai Jx- Jy \ 

A2 f 1 - f2 - Jz Jx + Jy Ai 

Al Jx + Jy ^2 — ^1 — Jz A2 

V Jx - Jy Al A2 -ei -62 + Jz J 



In order to facilitate our analysis of the decohering effect of the environment, in what follows we assume the 
Ohmic nature of the dissipativc reservoirs described by the spectral function S{uj) ~ J dte^'^'^ < ft^i, 2(0^1. 2(0) >= 
awcoth :^Q{u)c — w) of bandwidth uic (see the last part of this paper for a discussion of the relevance of this type of 
noise for superconducting qubits). 

The Ohmic nature of the noisy environment allows one to treat it in the standard Bloch-Rcdficld (i.e., weak-coupling 
and Markovian) approximation. In the basis spanned by the eigenstates of Eq.(3), the Bloch-Redfield equation for 
the reduced two-qubit density matrix reads [4]: 

Pnm{t) = —iijJnm Pnm{t) — Rnmkl Pkl{t), (4) 

k,l 

where cOnm = {En — Em) /Ji, n,m= 1, . . . , 4 are the transition frequencies between the eigenstates of Eq.(3), and the 
relaxation tensor 

Rnmkl = ^Im. ^ ^ ^nrrk + ^nk ^ ^ ^Irrm ~ ^Imnk ~ ^knml (^) 

r r 

is computed in terms of partial transition rates, whose real parts describe the effect of decoherence 

Re{Almnk} = -^S{Unk) [(^l,lm,'^l,nk + '^2,lm^^,nk] > (6) 

while Im{Kimnk} yield the Lamb shifts of the frequencies uJnk- 

As a customary measure of the coupled qubits' performance we use gate purity 
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where the density matrix p'{t) is evaluated for a complete set of 16 disentangled initial states: />'(0) = >< 4'^„|, 
= |*a >i 0|*6 >2,(a,6 1,- ••,4) which are composed of |*i >= | i>, |*2 >= | T>,l*3 (I i> +| T>)/V^, 
and 1*4 >= (I i> +i\ T>)/V2. 

Below, we restrict our analysis to the Hamiltonians that remain constant, H{t) = HQ6{t)9{to — t), for the entire 
duration to the gate operation, and, therefore, commute at any times {[H{t), H{t')] = 0). We then demonstrate 
that within this class of quasi-stationary Hamiltonians, the problem of constructing an optimized (coherence-wise) 
implementation of a given universal gate allows for a rather simple and physically transparent solution. 

To that end, we utilize the mechanism of decoherence suppression proposed in Ref. [5] where the decay rate of 
any initial state of an idling pair of coupled qubits ("quantum memory") was shown to be related to the spectral 
properties of the Hamiltonian. 

Namely, it was pointed out in Ref. [5] that the relaxation-related component of the overall decoherence can be 
significantly reduced by tuning the Hamiltonian parameters to the point in the parameter space where a pair of the 
lowest eigenvalues of the quasi-stationary Hamiltonian (3) becomes degenerate (see also the more recent Refs. [6] for 
corroborating results). 

The underlying (energy exchange-based) mechanism of the suppression of relaxation is based on the observation 
that near a degeneracy point and at low temperatures the partial relaxation rates (6) appear to be given by linear 
combinations of the transition frequencies J2i^j ^ijl^ijl where the coefficients Cij are essentially independent of Uij 
[5]. Therefore, the relaxation rates attain their minimum values at those points in the parameter space where the 
largest possible number of the transition frequencies vanish due to degeneracy. By contrast, a contribution from the 
other, pure dephasing, processes is generally unavoidable and can only be suppressed by lowering the temperature of 
the reservoirs. 

In order to further illustrate the above point, in Fig.l we plot the absolute value of the purity decay rate \dP/dt\ 
(which, when computed in the Markovian approximation, remains approximately constant at times shorter than the 
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decoherencc time) as a function of the components of the inter-qubit couphng Jy and Jz, while keeping the single-qubit 

fields Bi constant and the length of the vector J = + Jy + Jl fixed. On physical grounds, the latter constraint 

is justified by the fact that in any realistic qubit system an unlimited increase of J would eventually result in an 
unwanted leakage from the designated two-qubit Hilbert susbspace. 

The color plot of Fig.l demonstrates that \dPldt\ does attain its absolute minimum at the points characterized by 
the incidence of double degeneracy between the eigenvalues of Eq.(3) (see Eq.(15) below). Furthermore, albeit being 
less effective than its double counterpart, the onset of even a single degeneracy between the two lowest eigenvalues 
appears to provide a relative improvement, as compared to the generic (non-degenerate) situation. Even in this case, 
the two-dimensional degenerate ground state gets protected by the energy gap separating it from the rest of the 
spectrum that gives rise to the exponential suppression of relaxation at low temperatures. Considering the challenge 
of a high-precision tuning of the qubits' couplings in any realistic setup, the possibility of improving the gate quality 
with such a more relaxed constraint on the Hamiltonian parameters can be of a particular interest for practical 
realizations of robust quantum gates. 

Having identified the conditions required for a suppression of decoherence, we can now attempt to satisfy them in 
the cases of some popular universal gates. For a starter, we impose the less stringent condition of a single degeneracy 
between the lowest pair of energy levels, while looking for a one-step implementation of the standard CNOT gate 
^CNOT = diag{\\ 0][0 1 0] [0 1] [0 1 0]). 

In the noiseless case, this goal can be accomplished by directly computing the matrix-valued logarithm of the 
unitary matrix Xcnot- It turns out, however, that the latter features a high degree of ambiguity 



HcNOT = ilh In {Xcnot) = C{A + 3)0' 
Here to is the total duration of the gate operation 
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E is the unit matrix accounting for the possible "mismatch" phase 
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(observe that [A, B] = 0) and 
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while the matrices 
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(11) 



manifest the multi-valuedness of the matrix logarithm of the CNOT gate which can be obtained by exponentiating 
Eq.(lO) for an arbitrary choice of the integer (n^) as well as continuous (0, parameters in Eqs.(12) and (13), 
respectively. The existence of such an extensive invariant subspace within the equivalence class of CNOT (see below) 
is a special property of this particular gate. 

A straightforward analysis reveals that the CNOT gate can indeed be reproduced up to a global phase 0o = — 7r/4 
(so that Xcnot = cxp{—i7r/A — itoHo)) under the condition of a single spectral degeneracy {Ei = E2 = —0.875, E3 = 
0.625, E4 = 1.125) with the following choice of the parameters in Eq.(8): = = Jy = 0, A2 = 1.5, t\ = —0.25, £2 = 
Jz = —0.66 (hereafter all the energies are in units of 1/to). 

In Fig.2 we contrast the resulting optimized gate against the standard CNOT protocol consisting of a series of 
non-overlapping pulses (see, e.g., [7]) 



Xcnot = exp ( - 



^-^) exp ( - i^a^) 



exp ( - i-a^) exp ( - i-a^a^) exp ( - i- ^ J 



(12) 
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In order to compare the quality of the two protocols we solve the corresponding Bloch-Redficld equations (4) in the 
presence of weak (a << 1) noise in both cases. As a result, we find that the one-step implementation, where both 
the two- and one-qubit operations are carried out simultaneously, takes about 15% of the time of and appears to be 
approximately 10 times better (in terms of the gate purity) than the standard protocol (12). 

We mention, in passing, that one can also carry out a one-step implementation of other commonly used (e.g., RNOT 
and fast Fourier transform) gates by employing the Hamiltonians with single spectral degeneracies, thus significantly 
improving the gates' performances. As the next step, we attempt to impose the double degeneracy condition that, 
if at all attainable, is expected to result in an even better performance. It is, however, conceivable that this more 
stringent condition may not necessarily allow one to construct an arbitrary two-qubit gate. 

Nonetheless, even if this turns out to be the case, one can still elect to settle for a more readily achievable goal 
of constructing a gate X' which is equivalent (i.e. X' = U\ ® U^XU^ ® U2 where U\^2 and U[ 2 are single-qubit 
rotations), albeit not necessarily identical, to a target gate X. It should be noted, though, that in practice such a 
substitute solution might only be acceptable if all the one-qubit rotations can be performed sufficiently quickly, thus 
contributing negligibly towards the overall loss of coherence. 

Formally, the equivalence between a pair of gates is established on the basis of a coincidence of their Makhlin's 
invariants (see, e.g., [7]) 

^ ^ tr^[m{X)] ^ ^ tr^[m{X)]-tr[m\X)] 
^ lQDet[X\ ' ^ ADet[X] 

where m{X) = X^Xb, Xb = Q^XQ, and Q = diagi[l z][0 i 1 0][0 i - 1 0][1 - i]). 

A straightforward analysis shows that in the presence of a double spectral degeneracy Gi always turns out to be a 
real number, which fact severely restricts the set of gates that can be constructed this way. One example of a gate 
incompatible with the double degeneracy condition is provided by the VSWAP gate with Gi = — i/4. By contrast, 
the CNOT equivalence class (albeit not the original CNOT gate itself!) with the invariants Gi = 0, G2 = 1 can be 
readily attained under the condition of double degeneracy. 

Our theoretical findings can be of an immediate relevance to such a viable candidate to the role of robust two-qubit 
gates as a pair of the charge/flux Josephson junction qubits. Therefore, in what follows we specifically discuss this 
particular realization and focus on its most coherence- friendly regime (dubbed " quantronium" in Ref. [8]) where the 
individual qubits are tuned to their optimal points ei 2 = and coupled together both, capacitively [8] and inductively 
[9]. 

At the optimal point, the eigenvalues of Eq.(3) assume a particularly simple form 



Ei,2 = Jx T ^J {Ai + A2)^ + {Jy - J, 



^3,4 = - Jx ± ^(Ai - A2)2 + {Jy + J,)2 (14) 

Considering, for the sake of simplicity, the case of two identical qubits (Ai^2 = A), one can readily show that the 
double degeneracy condition {Ei = Ei,E2 = Es) is achieved at 

Jx = 0, JyJ. = A^ (15) 

Thus, in order to implement a transformation that belongs to the CNOT equivalence class under the double degeneracy 
condition it suffices to use a single rectangular interaction pulse with the parameters 

JCNOT = (0, ^iViJ^ + A2 - VJ' - A2), ^iViJ^ + A2 - VJ' - A2)) (16) 

From the Bloch-Redfield equations we find that corresponding protocol outperforms the CNOT-equivalcnt gate con- 
structed in Ref. [10] with the use of the Heisenberg inter-qubit coupling, as far as both the duration (« 2.3 shorter) 
and robustness (« 25 times lower purity decay rate) are concerned, thereby resulting in « 60 times lower overall 
purity loss 1 — P{ta)- 

For a large part, the uqibuitous CNOT gate owes its high popularity to the fact that a generic two-qubit gate 
requires no more than three CNOT applications complemented by local single-qubit rotations [11]. Recently, the 
authors of Ref. [7] put forward an alternative (dubbed as the "B") gate with the invariants Gi = G2 = 0. Unlike 
CNOT, however, the B-gate only needs to be used twice in order to implement an arbitrary two-qubit gate (up to 
single qubit rotations). 
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With an eye on the possibihty of a further optimization, we find that our approach which takes a full advantage 
of the double degeneracy condition offers a superior realization of the B-gate universality class, too. Namely, by 
choosing Ai^2 = I-, ei,2 = Jx = 0, Jy = 0.58, = 1.71 (in units of 1/to) we arrive at the one-step implementation of 
the equivalent of the B-gate that takes only « 56% of the duration of the protocol of Rcf. [7] and results in a much 
improved (f* 4.5 times slower) purity decay rate, thus decreasing the purity loss 1 — P{to) by a factor of « 8. 

Finally, we apply our theoretical analysis to the two-qubit Josephson junction setups of Ref. [12] with the parameters 
A ^ 10 GHz and J ^ 20 GHz. To that end, we first extract the pertinent value a ^ 0.01 of the dissipative coupling 
from the noise spectral function. The latter is simply proportional to the single-qubit relaxation rate, 1/Ti = ^S'(A), 
which was directly measured right at the optimal point and found to be of order 1/Ti w 0.1 GHz (see the last reference 
in Ref. [12]). 

With these data at hand, we estimate the purity decay which would occur in the course of implementing the above 
equivalent of the B-gate as 1 — Psito) ~ 0.03. However, this estimate would only hold for temperatures above ^ O.IK 
(albeit smaller than A), for at still lower temperatures dccohcrcncc was found to be dominated by the non-Ohmic 
(most likely 1//) noise which can not be treated the same way as the Ohmic one within the Bloch-Redfield approach. 

Correspondingly, in the case of the before mentioned equivalent of the CNOT gate one would obtain l—PcNOT{to) ~ 
0.15 which compares favorably to the result 1 — PcNorito) ~ 0.27 obtained in the second reference of Ref. [12]. 
However, our theoretical prediction still falls too short of the maximally acceptable error rate ~ lO"** which can be 
tackled with the existing error-correction codes. A further improvement would then require a coherence-conscience 
engineering of the environment, with the aim at significantly reducing the parameter a. 

It should be noted, however, that the authors of Ref. [12] attributed the loss of purity measured in that work solely 
to the errors in controlling the shape of the operation pulses (finite rise/fall time). For comparison, the sensitivity 
of our optimized protocols, in the presence of the infraction with reservoirs, with respect to variations 5xi of the 
control parameters appears to be reduced, thanks to the quadratic (as opposed to a generic linear in the case of not 
coherence- wise optimization) dependence upon such variations near the minima of the deviation functional (1). 

For example, we find that in order to keep 1 — P{to) under the 10"'' threshold the tolerance limit of (de)tuning of 
the optimized CNOT and B-gates' parameters has to be better than 0.3% of their values. 

Before concluding, it is worth mentioning that our approach differs from the previous proposals for constructing 
" supercoherent" qubits which, albeit being capable of providing an exponential suppression of decoherence, require 
at least four physical qubits which arc governed by the Hamiltonian that conserves total spin [13]. Apart from a 
larger number of physical qubits required to encode a logical one, the condition of spin rotational invariance of the 
Hamiltonian is unlikely to be fulfilled in any realistic solid-state qubits where (contrary to liquid-state NMR designs) 
the instantaneous values of the single- and two-qubit terms in Eq.(2) are often related, thus resulting in the local 
terms {B ^ 0) which break the spin-rotational invariance. 

To summarize, we present a systematic approach to constructing simplified (one-step) and optimized (highly robust) 
implementations of various two-qubit universal gates. The high stability against decoherence is achieved, thanks to 
the choice of the Hamiltonian parameters that provide for the degeneracies in the instantaneous energy spectrum, 
thereby resulting in the suppression of relaxation processes. Finally, we anticipate that the general nature of the 
decoherence suppression mechanism exploited in this work should make it possible to further generalize our results 
to the sequences of two-qubit as well as irreducible multi-qubit gates. 
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